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SUMMARY 


This  report  concerns  two  of  the  mechanisms  by  which  high- 
energy  protons  can  disable  silicon  semiconductors:  ionization 
and  displacement.  The  basic  difference  between  these  two 
mechanisms  is  fairly  concisely  explained  in  a  recent  text  as 
follows : 

"Displacement  effects  are  the  changes  that  occur  in  the 
properties  of  a  material  when  energetic  (here,  proton) 
radiation  displaces  (here,  silicon)  atoms  from  their 
normal  positions  in  the  lattice.  These  effects  are  in 
contrast  to  the  ionizaton  effects  that  are  the  result  of 
orbital  electrons  being  excited  away  from  their  parent 
(here,  silicon)  atoms  by  the  radiation"  (Ref.  1)  . 

The  reader  is  referred  to  that  text  for  a  detailed  discussion 
of  the  various  types  of  silicon  device  degradation  that  can 
result  from  proton  bombardment. 

Here  we  ask  the  question:  If  subjected  to  proton  bombard¬ 
ment,  would  a  silicon  semiconductor  device  first  fail  by 
excessive  ionization  or  by  excessive  d i splacement?  To  answer 
the  question,  we  follow  the  usual  custom  of  assuming  that  the 
two  disabling  mechanisms  act  independently  of  each  other. 

This  assumption,  in  turn,  permits  us  to  utilize  the  vast 
store  of  disablement  information  on  various  devices  that  has 
been  accrued  over  years  of  gamma-  and  neutron-irradiation 
experiments  (e.g.,  see  Ref.  2).  The  former  experiments 
(where  displacement  effects  are  negligible)  determine  total 


ionization-dose  failure  levels,  while  the  latter  (where  ioni¬ 
zation  effects  are  negligible)  determine  displacement-density 
failure  levels. 

To  estimate  ionization  doses  and  displacement  densities  all 
along  the  tracks  of  penetrating  protons,  one  of  course  needs 
at  least  a  reasonable  representation  of  the  proton  fluence  as 
a  function  of  both  the  position  and  energy  of  the  projectile 
particles.  By  following  the  type  of  approach  adopted  by 
Mather  and  Segre  in  1951  (Ref.  3),  such  a  fluence  function 
has  been  constructed  and  is  described  in  Reference  4.  This 
function  permits  an  immediate  estimate  of  the  ionization  dose 
along  the  path  followed  by  a  proton  beam  as  its  individual 
particles  straggle  to  their  various  stopping  locations  in  a 
thick  silicon  target. 

To  estimate  displacement  densities  along  the  path  one  also 
needs  a  second  function,  a  displacement  damage  coefficient 
which  is  dependent  on  proton  energy.  By  use  of  two  prelimi¬ 
nary  forms  of  this  second  function,  the  likely  magnitudes  of 
the  displacement  density  that  beams  of  25-  to  250-MeV  protons 
might  induce  in  thick  silicon  targets  have  been  explored 
(Ref.  5).  More  recently,  a  third  form  of  the  displacement 
damage  function  has  been  constructed.  As  fully  explained  in 
the  Appendix  of  this  report,  this  third  form  was  structured  to 
conform  with  a  wide  assortment  of  available  theoretical  and 
experimental  information.  In  the  following  sections,  we  will 
detail  the  findings  that  have  been  based  on  the  Appendix  dis¬ 
placement  damage  function  and  on  the  proton  fluence  mod(pl 
described  in  Reference  4. 


v-» 


Section  1  sets  forth  the  total  dose  profiles  in  thick  sili¬ 
con  targets  that  would  result  from  bombardments  by  beams  of 
25-  to  250-MeV  protons.  For  beams  of  an  incident  "reference 
fluence"  (10^3  protons/cm^,  corresponding  to  a  flow  of 
16  mA/m^  for  one  second) ,  the  dose  at  points  of  proton-entry 
ranges  from  ~l/2  Mrad  (250-MeV  protons)  to  ~3  Mrad  (25-MeV 
protons) .  The  peak  dose  (that  at  the  "Bragg  peak"  of  the 
dose  profile)  occurs  just  short  of  the  mean  penetration 
depth.  On  the  basis  of  the  Reference  4  fluence  model,  the 
peak  value  is  estimated  to  be  5.8  times  the  value  at  points 
of  proton-entry  into  the  thick  target.  This  factor  would  be 
appreciably  higher  if  straggling  were  not  taken  into 
account.  Where  the  dose  peaks,  the  fluence  of  surviving 
(i.e.,  unstopped)  protons  is  estimated  at  ~81%  of  the  inci¬ 
dent  fluence. 

In  Section  2,  we  estimate  the  capability  of  protons  (at  ad.1 
points  along  their  paths  in  the  thick  silicon  crystal  target) 
to  cause  silicon  atom  displacements  relative  to  the  same 
capability  of  1-MeV  neutrons.*  It  is  customary  to  relate  the 
proton- induced  displacement  density  to  that  induced  by  1-MeV 
neutrons,  since  the  displacement-disabling  fluences  of  such 
neutrons  have  been  measured  over  years  of  nuclear  weapons 
effects  experiments.  The  profiles  of  the  displacement  damage 
ratio  are  presented  for  beams  of  protons  with  incident 
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*Since  the  spectra  of  neutrons  from  sources  used  in  displace¬ 
ment  damage  experiments  vary  considerably,  it  has  become 
customary  to  equate  the  fluence  using  a  particular  source 
spectrum  to  the  "1-MeV-equivalent  neutron  fluence"  that  would 
cause  the  same  displacement  damage. 
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energies  in  the  same  25-  to  250-MeV  range  examined  in  our 
total  dose  study.  We  estimate  that  those  proton  beams  that 
are  capable  of  penetrating  at  least  ~10  g/cm^  of  target 
material  (i.e.,  with  incident  energies  >_  ~100  MeV)  would  be 
roughly  twice  as  effective  as  1-MeV  neutrons  in  causing  dis¬ 
placements  near  points  of  proton-entry.  The  peak  of  the 
displacement  density  profile  would  occur  in  the  near¬ 
neighborhood  of  the  mean  stopping  depth.  For  incident 
protons  of  100-MeV,  the  peak  proton-to-l-MeV-equivalent- 
neutron  displacement  damage  ratio  is  estimated  at  ~5 .  For 
250-MeV  incident  protons,  the  ratio  is  estimated  at  ~4 .  The 
peak  ratio  for  the  least  energetic  proton  beam  considered 
(25-MeV  incident  energy)  is  estimated  at  ~17 .  However,  the 
mean  penetration  depth  of  this  particular  beam's  protons 
would  be  only  ~0 . 8  g/cm^  of  target  material.  These  peak 
ratios  would  be  considerably  higher  if  straggling  did  not 
occur  . 

In  Section  3,  we  compare  the  irradiation  time  to  disable  a 
typical  bipolar  device  via  excessive  total  dose  to  that  for 
disablement  via  displacement  damage.  Reference  kill  criteria 
are  taken  at  1  Mrad(Si)  for  ionization  dose  disablement  and  at 
a  displacement  density  equal  to  that  resulting  from  exposure 
to  lO1-^  1-MeV  neutrons/cm^ .  Based  on  these  reference  kill 
criteria,  a  beam  that  delivers  a  flux  of  10^3  protons/cm^ *s 
to  a  thick  assembly  of  such  devices  could  always  first  kill 
each  such  device  by  ionization  alone.  Under  the  cited  con¬ 
ditions,  the  devices  located  near  the  mean  penetration  depth 
of  protons  would  fail  (via  ionization)  in  times  ranging  from 
”0.06  seconds  (25-MeV  beam)  to  0.3  seconds  (250-MeV  beam) . 
Devices  near  points  of  beam-entry  would  absorb  the  reference 
killing  dose  in  somewhat  longer  t  im.es:  0.3  seconds  for  the 

25-MeV  beam  and  2.0  seconds  for  the  250-MeV  beam .  Under 


reference  conditions,  the  time  required  to  cause  killing 
levels  of  the  displacement  density  is  estimated  to  be  a  fac¬ 
tor  of  ~4  to  ~12  larger  than  the  ionization  kill  time. 

We  should  here  remind  the  reader  that  the  ratio  of  ionization 
to-d i splacement  kill  times,  of  course,  scales  directly  with 
the  ionization  dose  kill  criterion  and  inversely  as  the 
displacement  kill  criterion.  For  example,  a  device  which 
fails  (via  displacements)  upon  exposure  to  only  10-*-^  1-MeV 
neutrons/cm^  or  to  our  reference  1-Mrad  ionization  dose, 
would  exper ience--essent ial ly  simultaneously--a  killing 
level  of  displacements  and  also  a  killing  level  of  ioniza¬ 
tion.  If  this  same  displacement-sensitive  device  was  so 
"total-dose-robust"  as  to  require  10  Mrads  for  ionization 
disablement,  it  would  first  be  killed  by  the  displacement 
mechanism  and,  considerably  later,  be  killed  again  by 
excessive  ionization. 

In  Section  4,  we  recall  our  original  plan  to  examine  the 
bombardment  of  silicon  semiconductors  by  deuterons  and 
tritons  as  well  as  by  the  proton  beams  considered  here.  The 
desirability  of  experiments  that  would  generate  the  type  of 
input  needed  to  extend  our  analyses  to  cover  all  three 
hydrogen- isotope  beams  is  noted.  Also,  our  plan  to  proceed 
analytically  with  similar  studies  involving  the  heavier 
isotope  beams  is  outlined. 

In  the  Appendix,  we  explain  the  construction  of  our  final 
form  of  the  displacement  damage  function. 
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SECTION  1 


PROTON  BEAM  TO  SILICON  TARGET 
ENERGY  TRANSFER  VIA  ATOM 
IONIZATION/EXCITATION  MECHANISMS 


Equation  (11)  of  Reference  4  permits  estimates  of  the  energy 
transferred  from  a  proton  beam  to  a  silicon  target  via  non¬ 
nuclear  interactions  (i.e.,  via  the  so-called  "dE/dx"  mecha¬ 
nisms  of  silicon  atom  ionization  and  excitation).  Reference  4 
considers  the  irradiation  of  a  uniform  semi-infinite  target 
by  a  uniform  monod i rect iona 1  beam  whose  axis  is  normal  to  the 
target  face.  The  energies  of  the  incoming  beam  particles  are 
highly  concentrated  about  an  energy,  Em,  which  is  one  of  the 
basic  input  quantities  that  define  a  beam/target  system.  A 
second  input  is  R,  the  mean  penetration  depth  of  beam  nuclei 
into  the  target  material.  For  the  proton/silicon  system  under 
study  here,  R  is  given  by 

R  =  3.05  x  10"3  x  Em  .  (1) 

Above,  b  =  1.74  and,  with  Em  in  MeV,  R  is  in  units  of  g/cm^ . 
The  third  input  is  r,  a  dimensionless  quantity  that  relates 
to  the  straggling  of  the  proton  stopping  points  in  the  thick 
silicon  target.  According  to  the  model  used  in  deriving 
Equation  (11)  of  Reference  4,  the  standard  deviation  of  the 
Gaussian  distribution  of  the  stopping  depths  is  R/r .  For 
protons,  the  model's  prescription  gives  r  =  71.4. 

Figure  1  shows  the  estimated  dE/dx  dose  that  would  be  felt  at 
points  of  a  thick  silicon  target  that  is  irradiated  by  beams 
that  deliver  a  reference  fluence  of  10i3  protons/cm^  to  the 
target  face.  To  accure  a  total  fluence  of  lO^-3  pa  r  t  i  c  1  es/cm^ 
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in  one  second,  the  average  current  density  at  the  target  must 
be  16  mA/m^. 

According  to  our  Figure  1  estimates,  a  reference  fluence  beam 
of  the  most  penetrating  particles  considered  (i.e.,  250-MeV 
protons)  would  deliver  ~l/2  Mrad  dose  to  silicon  near  the 
target  face  and  a  maximum  of  ~3  Mrad  at  the  Bragg  peak  of  the 
dE/dx  curve.  To  reach  silicon  at  depths  of  10  g/cm^  or  more, 
proton  energies  above  ~100  MeV  would  be  required.  And,  due 
to  the  energy  dependence  of  the  dE/dx  mechanisms,  it  is 
evident  that  the  lower  energy  beams  would  be  more  damaging  to 
the  silicon  devices  within  their  reach  than  would  be  the  more 
penetrating  higher  energy  beams  of  equal  fluence.  The  peak 
value  dose  is  a  factor  of  5.8  higher  than  the  near-surface 
dose  according  to  the  Reference  5  model  which  accounts  for 
straggling.  This  factor  would  be  significantly  higher  if 
straggl ing--which ,  of  course,  does  occur--did  not  occur. 
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SECTION  2 


SILICON  ATOM  DISPLACEMENT  DENSITIES: 
PROTON- INDUCED  RELATIVE  TO 
1-MeV- NEUTRON- INDUCED  DENSITIES 


The  energy  transferred  from  protons  to  silicon  atoms  via  the 
dE/dx  mechanisms  has  been  estimated  by  us  to  be  some  three 
orders  of  magnitude  larger  than  that  given  over  to  displacing 
these  atoms  from  their  normal  positions  in  the  silicon  crystal 
lattice.  However,  as  will  be  demonstrated  in  Section  3, 
dE/dx-  and  displacement-disablements  of  silicon  devices  may 
result  from  irradiations  that  are  within  only  one  order  of 
magnitude  of  each  other  in  duration.  To  lay  the  groundwork 
for  comparing  the  two  types  of  disablements  in  Section  3,  we 
here  first  relate  the  displacement  damage  capability  of  pro¬ 
tons  to  that  of  1-MeV  neutrons,  since  the  fluence  of  such 
neutrons  that  would  disable  various  semiconductors  via  the 
displacement  mechanism  is  known  (e.g.,  see  Ref.  2). 

As  customary  in  the  radiation  effects  field,  we  assume  that 
the  degradation  of  a  semiconductor's  performance  due  to 
displacements  is  a  function  of  the  displacement  number- 
density,  independent  of  the  type  of  radiation  that  causes  the 
displacements.  The  induced  displacement  density  is  taken  to 
be  proportional  to  the  density  of  primary  knock-on  atom  (PKA) 
displacements  times  an  average  "d i spl acemen t-e f f ec t i ve "  PKA 
recoil  energy,  w,  which  is  a  function  of  the  bombarding 
particle's  energy,  E.  As  fully  explained  in  the  Appendix,  the 
total  d i spl acemen t  density  (PKAs  plus  all  generations  of 
secondary  displacements)  at  depth  x  of  our  thick  silicon 
target  is  proportional  to 


(  J-  -)  •  ;(x,E)dE. 


Above,  i>(x,E)dE  represents  the  fluence  of  protons  at  depth  x 
with  energies  in  dE  about  E,  and  o  represents  the  displacement 
cross  section  for  these  protons.  Like  w,  o  is  a  function  of 
E.  The  <*>  ( x ,  E )  function  used  in  this  study  is  defined  by 
Equation  (8)  of  Reference  4  and  the  function  is  repre¬ 

sented  by  Curve  B  of  Figure  6  of  the  Appendix. 

The  curve  on  the  right  of  Figure  2  represents  that  portion  of 
the  complete  function  (the  above-cited  Curve  B)  which 

significantly  affects  our  displacement  estimates.  If  strag¬ 
gling  did  not  occur,  all  protons  would  travel  to  a  common 
depth  where  all  would  stop.  Also,  all  would  have  the  same 
energy  at  any  depth  along  the  way  to  the  common  stopping 
depth.  That  is,  the  displacement  density  at  any  depth  x 
would  simply  be  proportional  to  o(E)*w(E),  where  E  is  the 
common  energy  that  all  protons  would  have  at  this  depth. 

Since  the  function  increases  by  a  few  orders  of  magnitude 

as  proton  energy  is  reduced  by  dE/dx  mechanisms,  the  near- 
end-of-range  displacement  density  would  be  some  hundredfold 
larger  than  the  density  at  the  target  face,  j_f  straggling  did 
not  occur. 

To  help  understand  the  nature  of  the  displacement  density 
distribution  to  be  expected  when  straggling  is  taken  into 
account,  we  have  prepared  the  family  of  curves  shown  on  the 
left  of  Figure  2.  There,  the  solid-line  curves  give  the 
average  proton  energy  as  a  function  of  depth  in  the  silicon 
target,  this  average  being  estimated  by  use  of  Equation  (10) 
of  Reference  4.  The  dashed  curves,  obtained  by  use  of  Equa¬ 
tion  (9)  of  Reference  4,  indicate  the  percent  of  the  incident 
protons  that  remain  unstopped  at  each  depth.  For  the  displace¬ 
ment  density  to  be  high  at  a  given  depth,  we  would  expect  two 
conditions  would  nave  to  exist  at  that  depth.  First,  to  induce 


a  significant  number  of  displacements,  the  fluence  of  protons 
should  be  a  reasonable  percentage  of  the  incident  fluence. 
Second,  the  average  energy  of  the  protons  should  be  low  so 
that  the  quantity  o-u>r  to  which  the  displacement  density  is 
proportional,  would  be  high. 

Figure  2  permits  us  to  predict  the  general  shape  of  the 

displacement-density  profile  to  be  expected  when  a  beam  of 

protons  of  a  particular  incident  energy  bombards  a  thick 

silicon  target.  Starting  at  the  target  face  and  moving 

into  the  target  material,  we  note  that  this  density  should 

increase  initially  since  the  fluence  remains  essentially 

constant  and  the  average  proton  energy  drops  (i.e.,  o 

increases)  significantly.  As  we  move  to  depths  near  and 

beyond  the  Bragg  peak,  the  surviving  fluence  falls  off  so 

fast  as  to  overwhelm  the  increase  in  o*uj  that  results  from 

the  lessening  of  the  average  proton  energy.  Hence,  we  should 
« 

expect  that  the  displacement  density  would  go  to  a  maximum 
and  then  plunge  steeply  to  insignificant  levels  at  points 
slightly  beyond  the  Bragg  peak  neighborhood.  The  lower  the 
incident  proton  energy,  the  higher  the  expected  peak  of  the 
displacement  density,  since  lower  average  energies  will  be 
realized  by  a  significant  fraction  of  the  protons  in  this 
neighborhood.  However,  the  penetrability  of  the  lower  energy 
beams  is,  of  course,  very  low. 

In  Figure  3,  we  show  a  family  of  curves  that  have  profiles 
similar  to  the  displacement-density  profiles.  The  ordinate 
of  any  point  on  a  curve  represents  a  dimensionless  ratio. 

This  ratio  equals  the  value  of  the  integral  defined  in  the 
second  paragraph  of  this  section  divided  by  *0 • °N • ^  •  Here, 
f>0  is  a  fluence  equal  to  the  total  proton  fluence  at  the  face 
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of  our  thick  silicon  target;  a^(=  3.3  x  10”24  cm^)  represents 
the  scattering  cross  section  of  silicon  assumed  for  1-MeV 
neutrons  (Ref.  14);  and  0.6  x  [2  x  28/29^]  x  1  =  0.040 

MeV)  is  taken  as  the  average  recoil  energy  of  a  silicon 
nucleus  when  bombarded  by  1-MeV  neutrons.  The  0.6  factor  in 
the  bity  computation  is  that  cited  in  Reference  11  to  account 
for  the  forward  scattering  of  MeV  neutrons  in  the  center-of- 
mass  system.  These  assumed  data  give  =  1*33  x  10”25 

MeV^cm^  for  the  1-MeV  neutron  displacement  kerma. 

As  noted  in  Figure  3,  the  ordinate  compares  the  silicon-atom 
displacement  density  at  a  given  depth  per  proton/cm^  entering 
the  target  at  its  face  with  the  density  that  would  result  if 
the  silicon  at  that  same  depth  were  exposed  to  a  1-MeV  neutron 
fluence  of  1  neutron/cm^.  The  essential  point  here  is  that 
we  compare  displacement  densities  that  arise  from  a  uniform 
proton  bombardment  of  a  surface  (that  of  the  target  face) 
with  a  uniform  neutron  irradiation  of  a  volume  (that  of  the 
thick  silicon  target),  the  incoming  proton  fluence  and  the 
neutron  fluence  being  equal.  This  particular  ratio  readily 
lends  itself  to  the  accomplishment  of  our  primary  objective-- 
a  comparison  of  a  proton  beam's  capability  to  disable  a 
silicon  semiconductor  via  excessive  ionization/excitation 
with  its  capability  to  disable  via  excessive  displacements. 
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SECTION  3 


IONIZATION  VERSUS  DISPLACEMENT  DISABLEMENT 
OF  BIPOLAR  DEVICES 

According  to  the  state-of-the-art  hardness  levels  cited  in 
Reference  2,  typical  bipolar  technology  silicon  devices  can 
survive  an  average  total  ionizing  dose  of  1  Mrad  and  the 
effects  of  displacements  induced  by  exposure  to  an  average 
fluence  of  10^-^  1-MeV  neu  t  r  ons/cm^  .  For  reference  purposes, 
let  us  adopt  the  1-Mrad  figure  as  the  average  total  dose 
required  for  the  permanent  disablement  of  a  silicon  device. 
Also,  for  reference  purposes,  let  us  assume  that  the 
displacement-density  induced  by  exposing  a  device  to  an 
average  fluence  of  10^4  1-MeV  neutrons/cm^  is  sufficient  for 
permanent  disablement.  For  these  selected  reference  cri¬ 
teria,  we  now  compare  the  times  required  for  a  given  proton 
beam  to  kill  a  device  via  total  dose  with  that  required  for  a 
displacement  kill. 

In  Figure  4,  we  show  estimates  of  the  two  kill  times  for 
irradiations  by  a  beam  flux  of  10^  protons/cm^ -s .  The 
family  of  solid  curves  gives  the  estimated  total-dose  kill 
time  for  an  assembly  of  small  silicon  devices,  the  totality 
of  which  is  assumed  to  make  up  the  equivalent  of  a  thick 
silicon  target.  The  dashed  curves  give  the  estimated  times 
for  displacement  kills.  The  kill  time  for  a  device  that  is 
located  where  the  fluence  is  attenuated  to  0.1%  of  its  inci¬ 
dent  value  is  about  two  orders  of  magnitude  larger  than  that 
for  a  device  situated  at  any  depth  less  than  the  mean  proton 
ange.  This  is  a  reason  for  our  terminating  the  curves  at 
the  0.1%  fluence  points. 

As  seen  in  Figure  4,  which  is  based  on  the  reference  kill 
criteria,  disablement  via  total  dose  would  occur  prior  to 
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disablement  via  excessive  displacements.  In  Figure  5,  we 
show  the  ratio  of  the  dose-kill  time  to  the  displacement-kill 
time  for  the  reference  criteria.  We  conclude  that  kills  via 
total  dose  should  occur  in  ~1/10  to  ~l/5  of  the  times 
required  for  displacement  kills  for  devices  hardened  to  the 
levels  cited  on  Figures  4  and  5. 

It  may  be  noted  that  the  ratio  of  i on i za t i on- t o-d i splacemen t 
kill  times,  of  course,  scales  directly  with  the  ionization 
dose  kill  criterion  and  inversely  as  the  displacement  kill 
criterion.  For  example,  a  device  which  fails  (via  displace¬ 
ments)  upon  exposure  to  only  10^  1-MeV  neutrons/cm^  or  to  our 
reference  1-Mrad  ionization  dose,  would  exper i ence--essent i a  1 ly 
s imultaneously--a  killing  level  of  displacement  and  also  a 
killing  level  of  ionization.  If  this  same  displacement- 
sensitive  device  was  so  "total-dose-robust"  as  to  require  10 
Mrads  for  ionization  disablement,  it  would  first  be  killed  by 
the  displacement  mechanism  and,  considerably  later,  be  killed 
again  by  excessive  ionization. 
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SECTION  4 


CONCLUDING  REMARKS 

At  the  beginning  of  this  study,  we  planned  to  examine  the 
effects  of  silicon  bombardment  by  each  of  the  three  hydrogen 
isotopes  (i.e.,  by  deuterons  and  tritons  as  well  as  by  protons) 
Although  the  two  heavier  isotopes  cannot  penetrate  target 
material  as  deeply  as  can  protons  of  equal  energy,  they  are 
more  damaging.  For  example,  at  points  reachable  by  the 
heavier  isotopes,  both  the  dE/dx  energy  deposition  density 
and  the  displacement  number-density  will  be  higher  than  the 
corresponding  insults  induced  by  protons.  It  was  simple 
enough  to  estimate  the  dE/dx  energy  deposition  caused  by  each 
of  the  three  isotopes  by  use  of  the  Reference  4  fluence  model. 
The  total  dose  profiles  for  the  trio  of  projectiles  have,  in 
fact,  been  computed  on  the  basis  of  that  model  (Ref.  5). 
However,  we  have  thus  far  only  estimated  the  magnitudes  of 
peak  displacement  densities  to  be  expected  when  a  thick  sili¬ 
con  target  is  bombarded  by  each  of  the  three  isotopes  (Ref.  5). 

Now,  Rutherford  scattering  is  the  likely  dominant  displace¬ 
ment  mechanism  as  a  hydrogen  isotope  nears  the  end  of  its 
trip  into  a  thick  silicon  target.  And,  since  the  maximum 
displacement  density  occurs  near  the  journey's  end,  it  was 
possible  to  estimate  the  peak  displacement-density  for  each 
type  of  projectile  particle  by  use  of  the  Rutherford  model. 
However,  to  develop  the  complete  displacement  profiles,  one 
needs  a  reasonable  representation  of  the  displacement  damage 
coefficient  (i.e.,  the  function)  at  high  projectile 

energies  where  the  simple  Rutherford  model  is  not  applicable. 

In  the  case  of  proton  projectiles,  the  shape  of  the  o*.u  func¬ 
tion  over  much--but  not  all--of  the  high-energy  range  had 
been  established  by  the  Bell  Laboratory  proton  experiments 
(Ref.  6) .  As  explained  in  the  Appendix,  those  experiments-- 


together  with  the  material  presented  in  References  7,  8  and 
9--made  it  possible  for  us  to  construct  a  a*<*>  function  over 
the  full  high-energy  range  needed  for  the  proton- i nduced 
displacement  estimates  presented  here. 

At  this  writing,  we  do  not  have  experimental  data  for  the 
bombardment  of  silicon  by  high-energy  deuterons  or  tritons. 
Experiments,  similar  to  those  conducted  by  Bell  Laboratory 
workers  using  protons  within  a  limited  energy  range,  should 
be  performed  using  each  of  the  three  hydrogen  isotopes  over 
an  energy  range  reaching  to  ~300  MeV.  Such  thin-target 
experiments*  would  establish  the  shapes  of  the  functions 
for  each  of  the  three  projectile  types  in  the  complex  high- 
energy  regime,  a  regime  that — at  its  lower  end— over  laps  (and 
is  normalizable  to)  the  functions  in  the  Rutherford  scat¬ 
tering  regime. 

Pending  the  performance  of  the  above  types  of  experiments,  we 
plan  to  construct  functions  for  the  two  multiple-nucleon 

hydrogen  isotopes  (deuterons  and  tritons)  by  a  procedure  simi 
lar  to  that  used  for  the  simple  single-nucleon  projectile, 
the  proton.  That  procedure  (see  Appendix)  will  have  to  be 
modified  to  account  for  the  nearly-simultaneous  bombardment 
of  silicon  nuclei  by  two  or  three  nucleons,  those  initially 
assembled  in  the  rather  weakly  bound  configuration  that 
represents  a  deuteron  or  a  triton.  Armstrong  and  Colborn 
have  described  a  model  by  which  such  high-energy  bombardments 
may  be  treated  as  a  sequence  of  single  nucleon  interactions 
(Ref.  10) .  We  plan  to  examine  that  model  as  a  possible  means 
to  construct  the  o * ui  functions  for  deuteron  and  triton  beams. 

*  In  addition  to  thin-target  exposures  to  highly  monoener get i c 
projectiles,  thick-target  experiments  would,  of  course,  also 
be  desirable  since  these  could  reveal  the  smearing  effects  of 
straggling  by  direct  observation. 


With  these  functions  in  hand,  we  would  be  able  to  repeat  for 
the  heavier  hydrogen  isotopes  that  which  we  have  reported  here 
for  proton  beams.  This  would  complete  the  tasks  that  we  had 
in  mind  at  the  beginning--to  estimate  the  relative  vulnerabil¬ 
ities  of  silicon  semiconductor  to  disablements  via  excessive 
ionization  and  via  displacements  as  such  would  be  induced  by 
hydrogen  isotope  bombardment. 
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APPENDIX  : 


ASSUMED  REPRESENTATION  OF  THE 
DISPLACEMENT  DAMAGE  FUNCTION 

m 

Provided  that  a  proton-bombarded  silicon  atom's  recoil  energy 
is  not  too  high  but  high  enough  to  cause  its  displacement 
from  its  normal  position  in  the  crystalline  lattice,  the  total 
number  of  displacements  (primary  plus  all  secondaries)  has 
often  been  taken  to  be  proportional  to  the  recoil  energy  of 
the  primary  knock-on  atom  (e.g.,  see  p.  27  of  Ref.  11;  p.  104 
of  Ref.  12;  and,  p.  E-3  of  Ref.  13) .  However,  as  has  been 
noted  more  recently  in  References  7  and  8,  there  is  consider¬ 
able  experimental  evidence  that  the  energy  in  a  displacement 
cascade  does  not  increase  indefinitely  with  increasing  primary 
recoil  energy.  Thus,  for  the  recoil  energies  of  primary 
knock-on  atoms  (PKAs)  that  result  from  high-energy  proton 
collisions  with  silicon  nuclei  (i.e.,  where  proton  energy  Z 
coulomb  barrier  energy),  the  total  number  of  displacements 
will  be  proportional  to  an  average  "d  i  sp 1 acemen t -e f f ect i ve" 
recoil  energy,  which  is  less  than  the  average  PK A  energy. 

The  fluence  of  protons  with  energies  in  dE  about  E  at  depth 
x  in  a  silicon  target  is  denoted  by  i(x,E)dE.  This  incremen¬ 
tal  portion  of  the  total  fluence  at  x  will  cause  a  displace¬ 
ment  density  proportional  to  o  (  E)  •  -*>  (  E )  •  P(  x  ,  E )  d  E  .  Here, 

'i(E)  is  the  displacement  cross  section  for  protons  of  energy 
E,  and  *(E)  is  the  average  of  the  "d  i  sp  1  acemen  t- e  f  f  ect  i  ve" 
recoil  energies  for  such  protons.  We  refer  to  the  product 
■  *  *•  as  a  d  i  sp  1  acemen  t  damage  function  since  t  he  total  displace- 
m  ■  •  n  *■  -l-*ns  i  t  y  per  proton  is  taken  to  b>-  proportional  to  this 
*  a 'ic*  i  on  or  p  rot  on  enerjy.  As  n<<*  •  •  i  ubov  • ,  sines-  < 

1  *  n*  r  !  i  •  PK  A  r.-coi  1  ener  jy  1  ,  *  '•  "PKA  •  Figure  A 

i/s  ,;r  iS,siT--,j  r--pr>  n  t  it  ions  1  '  *  (f'-ir A;  i  nd 


{cmz  x  MeV) 


Basis  for  Curve  B 

(7)  Hard  sphere  silicon  atom 
model  joined  to  (2) 

@  Rutherford  model 

(3)  Bell  Lab  experimental 
displacement  data  fit  to  (2) 

(4)  Semi-empirical  fit  to  and 
extension  of  ©  (experimental 
cross-sections;  internuclear 
cascade  and  silicon 
displacement  cascade  machine 
calculations,  and,  approximate 
modeling  of  elastic  and 
nonelastic  proton-silicon 
collisions  I 


Proton  energv  (MeV) 

Frj'iru  6.  Assumed  representations  of  displacement  cross- 
section  1 imes  average  recoil  energy  of  primary 
knock-on  atoms  (Curve  A)  and  times  average 
"displacement-effective"  recoil  energy  (Curve  B) 
for  proton  bombardment  of  silicon  crystal. 
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As  indicated  on  the  figure,  our  o*w  function  was  constructed 
to  conform  with  available  experimental  data  plus  certain 
computed  information.  It  is  convenient  to  describe  the  con¬ 
struction  in  terms  of  four  bases,  denoted  by  (T),  (2),  (3),  and  (4) 
on  the  figure.  Basis  (3),  which  refers  to  the  19C4  experiments 
by  Bell  Laboratory  staff  (Ref.  6),  is  central  to  the  overall 
construction  of  Curve  B.  These  workers  irradiated  silicon 
solar  cells  by  protons  in  the  1-  to  130-MeV  range  and  deter¬ 
mined  the  dependence  of  the  displacement  damage  coefficient 
on  proton  energy  within  this  range.  They  observed  that--at 
the  beginning  of  the  range  i nves t iga ted-- the  coefficient 
decreased  with  increasing  proton  energy  as  predicted  by  the 
Rutherford  scattering  model.  Accordingly,  we  fit  the  exper¬ 
imental  Bell  curve  to  the  a,uipKA  curve  computed  by  use  of  the 
Rutherford  model  (see  overlapping  of  segments  (2)  and  (5)  i  n 
Figure  6 ) . 

According  to  the  prescription  of  Equation  2-8  of  Reference 
11,  we  assumed  the  Rutherford  model  to  apply  for  proton 
energies  down  to  ~1  KeV  and  the  hard  sphere  silicon  atom 
model  to  apply  for  energies  between  1  KeV  and  ~220  eV,  the 
range  of  curve  segment  (T).  The  220-eV  figure  represents  the 
lowest  energy  that  a  proton  can  have  and  still  displace  a 
silicon  atom  (30  eV  has  been  assumed  for  the  minimum  silicon 
recoil  energy  needed  for  a  displacement  to  occur).  The 
collision  cross  section  of  the  silicon  atom  in  the  segmentdT) 
energy  interval  is  assumed  to  vary  linearly  with  proton 
energy.  For  an  atom  at  rest,  we  assume  a  radius  equal  to 
the  Bohr  radius  of  the  atom's,  outer  electrons.  The  cross 
section  for  1-KeV  protons  is  adjusted  to  give  a  continuous 
J  *  -  function  where  segments  (T)  and  (2)  join. 

As  seen  in  Figure  2  of  the  text,  the  number  of  protons  that 
ire  slowed  to  an  average  energy  of  less  than  100  KeV  is 
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estimated  to  be  negligible  by  use  of  the  i>(x,E)  model  used  in 
this  study.  Hence,  our  estimates  of  displacement  densities 
cannot  be  sensitive  to  the  assumed  form  of  the  o*w  function 
in  and  near  the  segment-(l)  reg ime  of  Figure  6. 

Unfortunately,  the  130-MeV  upper  limit  of  the  proton-energies 
used  in  the  Bell  Laboratory  experiments  falls  short  of  the 
nominal  250-MeV  maximum  of  interest  to  us  here.  To  extend 
lurve  B  to  cover  the  desired  high-energy  range,  we  generated 
segment-(4}  which  runs  from  60  MeV  to  300  M^V  and  reasonably 
fits  the  Bell  data  in  the  60-  to  130-MeV  range.  For  proton 
energies  in  the  segment-(J)  range ,  °‘^pka  ^as  an  elastic  colli¬ 
sion  component,  OgM^e^PKA'  an<^  nonelastic  collison  component 
°n ‘  (Sn>  PKA • 

To  estimate  the  elastic  component,  we  assert  all  protons  are 
scattered  by  a  small  angle  equal  to  the  proton's  de  Broglie 
wavelength  divided  by  a  certain  diameter.  The  diameter  is 
that  of  a  circle  whose  area  equals  oQr  the  elastic  scattering 
cross  section  for  proton-silicon  collisions.  This  assertion 
y ields 

°e*  He)  PKA  =  11  •  (2) 

8M 

where  M  is  the  mass  of  the  silicon  nucleus  and  h  is  Plank's 
constant.  The  above  gives  0."?  x  10"24  Cm2 *MeV  for  the 
elastic  component. 

The  nonelastic  component  is  estimated  by  assuming  that  the 
proton  collides  elastically  with  one  of  the  nucleons  in  the 
silicon  nucleus  and  also  that  the  following  conditions  apply. 
Whichever  of  the  two  colliding  nucleons  emerges  from  the 
encounter  with  the  greatest  energy  is  assumed  to  move  on  a 


straight  line  and--without  further  encounters--to  leave  the 
nucleus.  The  nucleon  emerging  from  the  collision  with  the 
lesser  energy  is  assumed  to  be  stopped  within  the  nucleus. 

The  momentum  of  this  trapped  nucleon  goes  over  to  the  nucleus 
formed  by  the  entrapment.  The  energy  of  this  nucleon  goes  to 
( "n ) PKA '  to  excitation  energy,  E* ,  of  the  compound  system; 
and,  to  rotational  energy,  which  we  neglect  in  estimating  the 
nonelastic  component  of  °*^pkA‘  Where  these  conditions 
apply,  one  finds 

°n  *  (^n)PKA  =  °n  *  '  E*  ,  (3) 

M-m 

where  m  is  the  proton  mass  and  E*  is  the  average  value  of  E* . 
We  use  the  1968  Soviet  Monte  Carlo  computations  for  E*  [see 
o.  46  of  Ref.  9],  And,  for  lack  of  better  information  at 
this  time,  we  assume  the  nonelastic  cross  section,  °n ,  for 
proton  bombardment  of  silicon-28  equals  that  for  proton 
bombardment  of  aluminum-27.  Experimental  values  of  on  for 
aluminum  are  taken  from  page  9  of  the  same  Soviet  paper  just 
cited. 

The  sum  of  the  two  components  of  °,uJpka  thus  estimated  is 
represented  by  that  portion  of  Curve  A  in  Figure  6  that 
lies  in  the  proton  energy  interval  of  60  to  300  MeV.  To 
estimate  the  corresponding  portion  of  the  o*w  function  (i.e., 
Curve  B) ,  we  looked  to  References  7  and  8  for  guidance. 

As  previously  implied,  the  authors  of  those  two  papers  taught 
us  that--for  high-energy  recoils--the  "displacement-effective 
average  recoil  energy,  *>,  does  not  remain  proportional  to  the 
average  of  the  actual  PKA  recoil  energy,  1>PKA*  Rather,  the 


ratio  ",  -"pka  decreases  with  increasing  For  protons 

with  energies  between  60  and  130  MeV ,  we  could  fit  the  Bell 
Laboratory  experimental  findings  by  assuming  that  w  =  ^ka  if 
"PKA  ^  0.35  MeV  and  that  w  =  0.35  MeV  if  %>ka  ^  0.35  MeV. 

This,  together  with  the  assumption  that  silicon  nonelast ic--a 
well  as  ei  is  tic  —  cross  sections  equaled  those  for  aluminum 
(as  ['resented  on  p.  9  of  Ref.  9  )  ,  permitted  us  to  extend 
Curve  B  from  130-MeV  to  300-MeV  proton  energies  in  the  manner 
shown  in  Figure  6. 

With  the  0.35-MeV  limit  on  w ,  we  have  the  following  prescrip¬ 
tions  for  the  two  components  of  o*w  in  the  segment-(4)  energy 
interval. 

ae’^e  =  °*23  x  lO-24  cm2 • MeV ;  60  <  E  <  90  MeV,  (4a; 
oe*we  =  0.35  x  oe  cm2  *MeV ;  90  <  E  <  300  MeV,  (4b) 

and 

°n*'Ln  =  °*35  x  °n  cm2*MeV;  60  <  E  <  300  MeV.  (4c) 

Above,  E  represents  proton  energy  and  the  o’ s  are  implied  to 
be  in  units  of  cm2.  The  initial  fall-off  and  then  flattening 
of  Curve  B  for  E  >  60  MeV  result  from  the  similar  behaviors 
of  the  o's  as  functions  of  E  (see  p.  9  of  Ref.  9  )  . 
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ATTN:  YD 

SPACE  DIVISION/YE 
ATTN:  YE 

SPACE  DIVISION/YG 
ATTN.  YG 

SPACE  DIVISION/YN 
ATTN  YN 

STRATEGIC  AIR  COMMAND/INA 
ATTN:  INA 

STRATEGIC  AIR  COMMAND/NRI-STINFO 
ATTN:  NRI/STINFO 

STRATEGIC  AIR  COMMAND/XPFC 
ATTN:  XPFC 

TACTICAL  AIR  COMMAND/XPJ 
ATTN:  TAC/XPJ 

3416TH  TECHNICAL  TRAINING  SQUADRON  (ATC) 
ATTN.  TTV 

DEPARTMENT  OF  ENERGY 

DEPARTMENT  OF  ENERGY 
ATTN  ESHD 

UNIVERSITY  OF  CALIFORNIA 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN  L 13  D  MEEKER 
ATTN  L156JYEE 
ATTN  L  156  R  KALIBJIAN 
ATTN  L  53  TECH  INFO  DEPT  LIB 
ATTN  L 84 H KRUGER 
ATTN  WORVIS 
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DEPARTMENT  OF  ENERGY  (CONTINUED) 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN  D  LYNN 
ATTN  E  LEONARD 

SANDIA  NATIONAL  LABORATORIES 

ATTN  ORG  2 LOO  B  L  GREGORY 
ATTN  ORG  2126  J  E  GOVER 
ATTN  ORG  2144  PV  DRESSENDORFER 
ATTN.  ORG  2146  T  A  DELLIN 
ATTN  ORG  2150  J  A  HOOD 
ATTN  ORG  2320  J  H  RENKEN 
ATTN  ORG  2321  L  D  POSEY 
ATTN  ORG  5152  J  L  DUNCAN 
ATTN  TFWROBEL 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN  OSWR/NED 
ATTN.  OSWR/STD/MTB 

DEPARTMENT  OF  TRANSPORTATION 


ATTN 

ARD  350 

NASA 

ATTN: 

CODE  313  V  DANCHENKO 

ATTN 

CODE  600  E  STASSINOPOULOS 

ATTN. 

CODE  660  J  TRAINOR 

ATTN 

CODE  695  M  ACUNA 

ATTN 

CODE  724  1  M  JHABVALA 

NASA 

ATTN: 

M  BADDOUR 

NASA  HEADQUARTERS 

ATTN 

CODE  DP  B  BERNSTEIN 

NATIONAL  BUREAU  OF  STANDARDS 

ATTN 

C  WILSON 

ATTN 

CODE  A327  H  SCHAFFT 

ATTN 

CODE  A347  J  MAYO  WELLS 

ATTN 

CODE  A353  S  CHAPPELL 

ATTN 

CODE  A361  J  FRENCH 

ATTN 

CODE  C216  J  HUMPHREYS 

ATTN 

T  RUSSELL 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ADVANCED  RESEARCH  &  APPLICATIONS  CORP 

ATTN 

L  PALKUTI 

ATTN 

R  ARMISTEAD 

AEROJET  ELEC 

TRO  SYSTEMS  CO 

ATTN 

D  TOOMB  B  1  Ml  D  4  34  1 

ATTN 

P  ..ATHROP 

AEROSPACE  CORP 


ATTN 

C  E  BARNES 

ATTN 

D  FRESH 

ATTN 

D  SCHMUNK 

ATTN 

G  GILLEY 

ATTN 

1  GARFUNKEL 

ATTN 

J  REINHEIMER 

ATTN 

J  STOLL 

ATTN 

J  WIESNER 

ATTN 

J  B  BLAKE 

ATTN 

K  T  WILSON 

ATTN 

M  DAUGHERTY 

ATTN 

N  SRAMEK 

ATTN 

P  BUCHMAN 

ATTN 

R  SLAUGHTER 

ATTN 

W  CRANE 

ATTN 

W  KOLASINSKI 

ALLIED  CORP 

ATTN 

DOCUMENT  CONTROL 

ALLIED  CORP.  BENDIX  FLIGHT 
ATTN:  E  MEEDER 

AMPEX  CORP 

ATTN:  K  WRIGHT 
ATTN  PPEYROT 

ANALYTIC  SERVICES.  INC  (ANSER) 

ATTN  ASHOSTAK 
ATTN  J  OSULLIVAN 
ATTN  P  SZYMANSKI 

APPLIED  SYSTEMS  ENGRG  DIRECTOR 

ATTN  J  P  RETZLER  NUC  S/V  MANG 

AVCO  SYSTEMS  DIVISION 
ATTN  D SHRADER 

BDM  CORP 

ATTN  C  M  STICKLEY 

BDM  CORP 

ATTN  DWUNSCH 

BEERS  ASSOCIATES.  INC 
ATTN  B BEERS 

BOEING  CO 

ATTN  M  ANAYA 
ATTN  C ROSENBERG 
ATTN  A  JOHNSTON 
ATTN  DEGELKROUT 
ATTN  ELSMITH 
ATTN  I  ARIMURA 
ATTN  R CALDWELL 
ATTN  W  DOHERTY 
ATTN  H  WICKLEIN 


DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 


BOEING  CO  (CONTINUED) 

ATTN  OMULKEY 
ATTN  C  DIXON 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 
ATTN  W  PRICE 

CAlSPAN  CORP 

ATTN  R  THOMPSON 

CHARLES  STARK  DRAPER  LAB.  INC 
ATTN  J  BOYLE 
ATTN  N  TIBBETTS 
ATTN  PGREIFF 
ATTN:  W  D  CALLENDER 

CINCINNATI  ELECTRONICS  CORP 
ATTN  L  HAMMOND 

CLARKSON  COLLEGE  OF  TECHNOLOGY 
ATTN  PJ  MCNULTY 

COMPUTER  SCIENCES  CORP 
ATTN  ASCHIFF 

Denver,  university  of 

ATTN  SEC  OFCR  FOR  F  VENDITTI 

DEVELCO.  INC 

ATTN  G  HOFFMAN 

E  SYSTEMS.  INC 

ATTN  K  REIS 

E  SYSTEMS.  INC 

ATTN  DIV  LIBRARY 

EATON  CORP 

ATTN  R  BRYANT 

ELECTRONIC  INDUSTRIES  ASSOCIATION 
ATTN  JKINN 

EMC  CORP 

ATTN  L STATES 

FORD  AEROSPACE  &  COMMUNICATIONS  CORP 
ATTN  TECH  INFOR SVCS 

GARRETT  CORP 

ATTN  S  ELLIOTT 

GENERAL  ELECTRIC  CO 

ATTN  DEDELMAN 
ATTN  DTASCA 
ATTN  DOCUMENTS  LIB 
ATTN  HODONNELL 


ATTN:  J  ANDREWS 
ATTN:  R  BENEDICT 
ATTN:  R  CASEY 
ATTN:  TECHNICAL  LIB 

GENERAL  ELECTRIC  CO 

ATTN:  B  FLAHERTY 
ATTN:  G  BENDER 
ATTN:  L  HAUGE 

GENERAL  ELECTRIC  CO 

ATTN:  G  GAT  I  MD 

GENERAL  ELECTRIC  CO 

ATTN:  C  HEWISON 
ATTN.  D  COLE 

GENERAL  ELECTRIC  CO 
ATTN:  J  MILLER 

GENERAL  RESEARCH  CORP 
ATTN.  A  HUNT 

GftORGE  WASHINGTON  UNIVERSITY 
ATTN:  A  FRIEDMAN 

GOODYEAR  AEROSPACE  CORP 

ATTN:  SECURITY  CONTROL  STA 

GRUMMAN  AEROSPACE  CORP 
ATTN:  J ROGERS 

GTE  GOVERNMENT  SYSTEMS  CORPORATION 
ATTN:  J A  WALDRON 

HARRIS  CORP 

ATTN:  E  YOST 
ATTN:  WABARE 

HARRIS  CORP 

ATTN  J  W  SWONGER 

HONEYWELL,  INC 

ATTN:  DHEROLD 
ATTN  D  LAMB 
ATTN  D  NIELSEN 
ATTN  R  BELT 
ATTN  RGUMM 

HONEYWELL.  INC 

ATTN  J SCHAFER 
ATTN  MS  7?5 5 

HUGHES  AIRCRAFT  CO 

ATTN  RC  HENDERSON 

HUGHES  AIRCRAFT  CO 

ATTN  W  SCHENET 


29 


DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 

HUGHES  AIRCRAFT  CO 
ATTN  J  HALL 

HUGHES  AIRCRAFT  COMPANY 
ATTN  A NAREVSKY 
ATTN  E  KUBO 
ATTN  L  DARDA 

IBM  CORP 

ATTN  H  MATHERS 

IBM  CORP 

ATTN  J  ZIEGLER 

IBM  CORP 

ATTN  AEDENFELD 
ATTN  N  HADDAD 

1 1 T  RESEARCH  INSTITUTE 
ATTN  AKBUTI 
ATTN  I  MINDEL 

ILLINOIS  COMPUTER  RESEARCH.  INC 
ATTN  ES  DAVIDSON 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN  TECH  INFO  SVCS 

IRT  CORP 

ATTN:  JAZAREWICZ 
ATTN  J  HARRITY 
ATTN.  J  M  WILKENFELD 
ATTN  M  ROSE 
ATTN:  MDC 
ATTN  N  RUDIE 
ATTN  R  MERTZ 

JAYCOR 

ATTN.  DR  M  TREADAWAY 
ATTN  R STAHL 
ATTN  T  FLANAGAN 

JAYCOR 

ATTN  R  SULLIVAN 

JAYCOR 

ATTN  C  ROGERS 
ATTN  R  POLL 

JOHNS  HOPKINS  UNIVERSITY 
ATTN  P  PARTRIDGE 
ATTN  R  MAURER 

JOHNS  HODKINS  UNIVERSITY 
ATTN  G  MASSON 


KAMAN  SCIENCES  CORP 
ATTN  C  BAKER 
ATTN  DIR  SCI  &  TECH  DIV 
ATTN  JERSKINE 
ATTN  N  BEAUCHAMP 
ATTN  W  RICH 

KAMAN  SCIENCES  CORP 
ATTN  E CONRAD 

KAMAN  SCIENCES  CORPORATION 

ATTN  TECH  LIB  FOR  D  PIRIO 

KAMAN  TEMPO 

ATTN  DASIAC 
ATTN  R  RUTHERFORD 
ATTN:  W  MCNAMARA 

KAMAN  TEMPO 

ATTN:  DASIAC 

LITTON  SYSTEMS,  INC 

ATTN:  EL  ZIMMERMAN 
ATTN:  F  MOTTER 

LOCKHEED  MISSILES  &  SPACE  CO.  INC 
ATTN  FJUNGA 
ATTN:  J  SMITH 
ATTN:  REPORTS  LIB 

LOCKHEED  MISSILES  &  SPACE  CO.  INC 
ATTN:  BKIMURA 
ATTN:  EHESSEE 
ATTN:  J  C  LEE 
ATTN  JCAYOT 
ATTN:  L  ROSSI 
ATTN.  P  BENE 
ATTN:  S  TAIMUTY 

LTV  AEROSPACE  &  DEFENSE  COMPANY 
ATTN  A  R  TOMME 
ATTN  LIBRARY 
ATTN:  TECH  DATA  CTR 

MAGNAVOX  ADVANCED  PRODUCTS  &  SYS  CO 
ATTN:  W  HAGEMEIER 

MARTIN  MARIETTA  CORP 
ATTN  J  TANKE 
ATTN.  J  WARD 
ATTN  MP  1601  W  BRUCE 
ATTN  RGAYNOR 
ATTN  TIC/MP  30 

MARTIN  MARIETTA  CORP 
ATTN  T  DAVIS 


>.y. 
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DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 


MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN:  M  POLZELLA 
R  ANDERSON 
R  KASE 

RESEARCH  LIB 


NORDEN  SYSTEMS.  INC 

ATTN  N  RIEDERMAN 
ATTN  TECH  LIBRARY 


ATTN 

ATTN 

ATTN 


NORTHROP CORP 

ATTN  A  BAHRAMAN 
ATTN  J  SROUR 
ATTN  ZSHANFIELD 


MARYLAND.  UNIVERSITY  OF 
ATTN.  H  C  LIN 


MCDONNELL  DOUGLAS  CORP 


NORTHROP CORP 

ATTN:  E  KING 
ATTN:  S  STEWART 


ATTN 

ATTN 

ATTN 

ATTN 


A  P  MUNIE 
D  L  DOHM 
M  STITCH 
R  L  KLOSTER 


PACIFIC  SIERRA  RESEARCH  CORP 

ATTN  H  BRODE.  CHAIRMAN  SAGE 


MCDONNELL  DOUGLAS  CORP 
ATTN  P ALBRECHT 


PHYSICS  INTERNATIONAL  CO 
ATTN:  J  SHEA 


MCDONNELL  DOUGLAS  CORP 
ATTN:  TECH  LIB 


POWER  CONVERSION  TECHNOLOGY.  INC 
ATTN:  V  FARGO 


R  &  D  ASSOCIATES 


MESSENGER.  GEORGE  C 

ATTN 

B  LAMB 

ATTN:  G  MESSENGER 

2  CYS  ATTN 

E  MARTINELLI 

2  CYS  ATTN 

G  MESSENGER 

MISSION  RESEARCH  CORP 

2  CYS  ATTN 

G  SAFONOV 

ATTN:  C  LONGMIRE 

ATTN 

M  GROVER 

MISSION  RESEARCH  CORP 

ATTN 

W  KARZAS 

ATTN:  R  PEASE 

RAND  CORP 

ATTN 

C  CRAIN 

MISSION  RESEARCH  CORP 

ATTN:  J  LUBELL 

ATTN 

P  DAVIS 

ATTN:  R CURRY 

RAND  CORP 

ATTN:  W  WARE 

ATTN 

B  BENNETT 

MISSION  RESEARCH  CORP.  SAN  DIEGO 

RAYTHEON  CO 

ATTN:  J  RAYMOND 

ATTN 

G  JOSHI 

ATTN:  WAN  LINT 

ATTN 

J  CICCIO 

MITRE  CORPORATION 

RAYTHEON  CO 

ATTN:  M  FITZGERALD 

ATTN 

A  VAN  DOREN 

MOTOROLA.  INC 

ATTN 

H  FLESCHER 

ATTN  A  CHRISTENSEN 

RCA  CORP 

ATTN 

G  BRUCKER 

MOTOROLA.  INC 

ATTN:  C  LUND 

ATTN 

V  MANCINO 

ATTN  L CLARK 

RCA  CORP 

ATTN:  0  EDWARDS 

ATTN 

R  SMELTZER 

NATIONAL  SEMICONDUCTOR  CORP 

RCA  CORP 

ATTN  FC  JONES 

ATTN 

E  SCHMITT 

NEW  MEXICO.  UNIVERSITY  OF 

ATTN 

W  ALLEN 

ATTN  DR NEAMEN 

RCA  CORP 

ATTN 

E  VAN  KEUREN 

31 


f  V, 


m 


& 


’>1 


m 


o-j 


i 


o 


lOVJ 


I  j  w 


:>u 


•.tv 


a 


itS 


DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 

RENSSELAER  POLYTECHNIC  INSTITUTE 
ATTN  RGUTMANN 

RESEARCH  TRIANGLE  INSTITUTE 
ATTN  M  SIMONS 

ROCKWELL  INTERNATIONAL  CORP 
ATTN  A  ROVELL 
ATTN  GA50  TIC/'L  G  GREEN 
ATTN  J  BELL 
ATTN  J  BURSON 
ATTN  VOE  MARTINO 
ATTN  V  STRAHAN 

ROCKWELL  INTERNATIONAL  CORP 
ATTN  TICAJOl 

ROCKWELL  INTERNATIONAL  CORP 
ATTN  LW  PINKSTON 
ATTN  TIC  124-203 

ROCKWELL  INTERNATIONAL  CORP 
ATTN  T  YATES 
ATTN  TICBA08 

ROCKWELL  INTERNATIONAL  CORP 
ATTN  J  BURSON 

SCIENCE  APPLICATIONS  INTI  CORP 
ATTN  0  LONG 
ATTN  DMILLWARD 
ATTN  D  STROBEL 
ATTN  L SCOTT 
ATTN  R  FITZWILSON 
ATTN  RJBEYSTER 
ATTN  V  ORPHAN 
ATTN  V  VERBINSKI 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN  J SPRATT 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN  W  CHADSEY 

SCIENCE  APPLICATIONS  INTL  COPP 
ATTN  PAZIELIE 

SCIENTIFIC  RESEARCH  ASSOC.  INC 
ATTN  HGRUBIN 

SINGER  CO 

ATTN  RSPIENGEL 
ATTN  TECH  INFO  CNTR 

SPERRY  CORP 

ATTN  J  INDA 


SPERRY  CORP 

ATTN:  P  MARROFFINO 

SRI  INTERNATIONAL 

ATTN:  A  PADGETT 

SUNDSTRAND  CORP 

ATTN  C  WHITE 

SYSTEM  DEVELOPMENT  CORP 

ATTN.  PRODUCT  EVAL  LAB 

SYSTRON  DONNERCORP 
ATTN:  J  RAY 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  G  R  EZELL 

TELEDYNE  SYSTEMS  CO 
ATTN:  RSUHRKE 

TEXAS  INSTRUMENTS,  INC 
ATTN:  E  JEFFREY 
ATTN:  F  POBLENZ 
ATTN:  J  SALZMAN 
ATTN:  T CHEEK 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 


ATTN 

A  WITTELES 

ATTN 

D CLEMENT 

ATTN 

F  FRIEDT 

ATTN 

H  HOLLOWAY 

ATTN 

M  S  ASH 

2  CYS  ATTN 

0  ADAMS 

ATTN 

PGUILFOYLE 

ATTN 

PR  REID 

2 CYS  ATTN 

R  PLT  BUCH,  HARD  &  SURV  LAB 

ATTN 

R  VON  HATTEN 

ATTN 

TECH  INFO  CTR.  DOC  ACQ 

ATTN 

W  B  ADELMAN 

ATTN 

W  ROWAN 

ATTN 

W  WILLIS 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  C  BLASNEK 
ATTN:  F  FAY 
ATTN  J GORMAN 

VISIDYNE.  INC 

ATTN  CH  HUMPHREY 
ATTN  W  P  REIDY 

WESTINGHOUSE  ELECTRIC  CORP 
ATTN  D  GRIMES 
ATTN  H  KALAPACA 
ATTN  R  C.RICCHI 

WESTINGHOUSE  El  ECTRIC  CORP 
ATTN  SWOOD 


